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We investigate theoretical approaches to pion-nucleus elastic scattering at high 
energies (300 < T n < 1 GeV). A "model-exact" calculation of the lowest-order 
microscopic optical model, carried out in momentum space and including the full 
Fermi averaging integration, a realistic off-shell pion-nucleon scattering amplitude 
and fully covariant kinematics, is used to calibrate a much simpler theory. The sim- 
pler theory utilizes a local optical potential with an eikonal propagator and includes 
the Coulomb interaction and the first Wallace correction, both of which are found to 
be important. Comparisons of differential cross sections out to beyond the second 
minimum are made for light and heavy nuclei. Particularly for nuclei as heavy as 
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Ca, the eikonal theory is found to be an excellent approximation to the full theory. 
PACS number(s): 25.80.Ek, 24.10.Eq, 24.10.Jv 
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I. INTRODUCTION 



The pion-nucleus interaction on and below the A 33 resonance has been extensively stud- 
ied, and microscopic models of elastic scattering do a reasonable job JE] of describing the 
data. The pion-nucleus interaction at energies above the A 33 resonance has received 0,|3],[4|,|5j 
much less attention. In this energy region (300 MeV < < 1 GeV), the pion has a much 
shorter wavelength. For example, the wavelength at resonance is about 4 fm (about the size 
of the nucleus) while at 1 GeV the wavelength is 1 fm (about the size of a single nucleon). 
The shorter wavelength implies that elastic and inelastic data at the higher energies will pro- 
vide sensitivity to the details of the spatial dependence of the ground-state and transition 
densities. The data may also prove sensitive to modifications |||7|] of the properties of the 
nucleon in the nuclear medium, or reveal significant contributions from exchange currents 

Moreover, the pion-nucleon two-body interaction becomes much weaker as one goes to 
energies above the A 33 resonance. The two-body total cross section becomes less than 30 
mb, which is about 15 % of that on the A 33 resonance. The first implication of this weaker 
amplitude is that the pion is able to penetrate deeper into the nucleus. A simple estimate 
H from total reaction cross section studies shows that a projectile can penetrate into a 
target to a radius that is equal to the impact parameter at which the profile function equals 
one mean free path. In Fig. [I] the profile function for 12 C, 40 Ca and 208 Pb are pictured, 
and the arrows indicate approximately how far into a nucleus a pion of the labeled energy 
can penetrate. The pion in the energy region from 500 MeV to 1 GeV is one of the most 
penetrating of the strongly interacting particles. 

The second implication of the weaker two-body cross section is that multiple scattering 
theory for the optical potential becomes increasingly convergent. A simple estimate of the 
convergence is obtained by comparing a typical second-order term in the optical potential 
to the first-order term. For the case of short-range correlations, one obtains |10| 



where a is the total two-body cross section, i c is the correlation length, and p is the nuclear 
density at which the pion interacts (See Fig. 1). The factor of 1/k, where k is the incident 
pion momentum, reflects the suppression of the pion propagator in with increasing 
energy. On resonance, we find R ~ 0.1, with R < 1 only because p is so small there. 
At 500 MeV, we find R « 0.04 and at 1 GeV, R « 0.02. Thus differences between the 
data and a carefully calculated result utilizing a first-order optical potential would be a 
strong indication of the presence of unconventional phenomena such as mesonic current 
contributions, modified nucleon properties in the medium, or other yet to be thought of 
effects. 

One option for calculating high-energy pion scattering is a momentum-space optical- 



potential approach. This presents the opportunity [11] to calculate the scattering from a 



lowest-order optical potential including the following features 1.) exact Fermi-averaging 



integration, 2.) fully covariant kinematics |12[], normalizations, and phase-space factors, 3.) 



invariant amplitudes [0,0 and 4.) finite-range, physically motivated two-body off-shell 
jn| amplitudes. Within the multiple scattering theory developed in Refs. [|I|,|T4] this is an 
exact calculation of the lowest-order optical potential. A brief review of this approach will 
be given in Sec. |TJ. 

A practical difficulty arises when one wants to use the momentum-space optical model to 
study high-energy pion-nucleus scattering. As the projectile energy increases, more partial 
waves are needed in the pion-nucleon two-body amplitudes. Below 300 MeV, only S-waves 
and P- waves are needed. Above 300 MeV, D-waves become important. F-waves become 
significant above 500 MeV and G-waves and H-waves above 700 MeV. At the same time 
the number of pion-nucleus partial waves is increasing at a rate proportional to the pion 
momentum. At the high energies the momentum-space approach becomes prohibitively 
computer-intensive, and one would like to search for a simpler alternative. 

The semi-classical theory immediately comes to mind as an alternative to the 
momentum-space approach. It is often used at high energies not only because it is eas- 
ier to compute, but also because the simpler character of the theory facilitates obtaining 



physical insights into the reaction. However, the semi-classical theory is expected to be a 
good approximation only for local potentials and only when the wavelength is sufficiently 
short. Because the exact pion-nucleus optical potential appears to be highly nonlocal, we 
examine the role that nonlocalities play in high-energy pion-nucleus scattering. This is done 
numerically in momentum space and is presented in Sec. |11 A| . We undertake this study as 
a first step in obtaining a more quantitative measure of the validity of the semi-classical 
theory than currently exist in the literature. 

In Sec. |l IB] we present our semi-classical model |L5]. In order to have a quantitative 
as well as a simple model for confronting data, we include the Coulomb interaction and 



the first Wallace correction |16 |. We establish the validity of the semi-classical model both 
by examining the size of the Wallace correction and by comparing it to the model-exact 
results obtained in momentum space. Comparisons of the eikonal model to model-exact 



calculations are made in Sec. Ill 



A previous similar investigation can be found in Ref . ]T7[ . There a factorized approxima- 
tion was used in the momentum-space calculations, so that the full nonlocality of the 
optical potential was not considered. Furthermore, Glauber multiple scattering theory JT9| 
was proposed, whereas we will examine an even simpler model, that of a local optical po- 
tential with the scattering solved via an eikonal approximation. A more detailed discussion 



that compares and contrasts our results with previous work is given in Sec. [TV]. A summary, 
conclusions and future prospects are presented in Sec. [V|. 



II. SCATTERING THEORY 

At the higher-energies where the pion-nucleon amplitude becomes a smooth function of 
energy and is much weaker than on resonance, the first-order optical potential may well be 
adequate for describing all of the conventional nuclear physics phenomena that enter into the 
reaction dynamics. Whether this is so remains for further investigation. Independent of this 
question, there is the question of whether simpler approaches to the calculation of the lowest- 
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order optical potential and the scattering from this potential are reasonably quantitative. 
It is this question upon which we concentrate here. Below we first review the momentum- 
space optical potential which serves as our 'model-exact' calculation. In momentum-space, 
we examine the importance of the nonlocalities caused by both the resonance propagation 
and the finite range of the two-body amplitude. Finding that these both can be accurately 
approximated, we present a simple eikonal model. 

Both the optical and the eikonal models use the same target wave functions, which 



are obtained from Hartree-Fock calculations p0|j21[| with spurious center-of-mass motion 
removed as discussed in A. They also use the same on-shell pion-nucleon two-body ampli- 
tudes, which are constructed from Arndt's [22] and Holder's p3| phase-shifts. 



A. Optical potential 

A complete description of the momentum-space optical potential that we use can be 
found in Ref. ITD, and a detailed discussion of the covariant kinematics that we use can be 



found in Refs. [12|]. The formal multiple-scattering theory in which this work is embedded 
is given in Refs. 0,0. Here we provide only a brief overview. 

The first-order optical potential in the impulse approximation can be written as, 



$A i uue) I U.) = £ / f£i |f f£ (*L, S , I U^) 

x {k% | t(E) | k n k N ) {k N k A -i | V^J , (2) 



where k^, k^ and k A -i are the momenta of the pion, the struck nucleon and the A — 1 
residual nucleons in the pion-nucleus center-of-momentum frame respectively. E is the 
incident energy of the pion and the nucleus (in the pion-nucleus center-of-momentum frame 
also), and a is a set of quantum numbers that specifically label the nuclear bound state. 



The target wave function (k^k^-i | ^ a £ ) defined covariantly [IT] contains a momentum- 
conserving delta function 5{kA — k^ — fc^-i) and is a function of the relative momentum 
between the nucleon and the A — 1 residual nucleons. The pion-nucleon t-matrix similarly 
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contains a momentum-conserving delta function and is a function of the relative momentum 
between the pion and the nucleon, k, 

(k' N k' n | t(E) | k N K) = 5(k' N + k' w -k N - K)(k' I t[u a (E, k n , k N )} | k) . (3) 

One of the these three momentum-conserving delta functions leads to overall momentum 
conservation; the others allow one to perform two of the three integrals in Eq. |2|. The 
remaining integral is the Fermi averaging integration which must be performed numerically. 
The details of how we perform the integration can be found in ||1 1|| . 

The energy at which one evaluates the two-body t-matrix, u> a (E, k^, k^), must be care- 



fully chosen |24| if a convergent perturbation theory is to result for calculations at and below 
the A33 resonance region. The energy u a is defined covariantly by first defining the energy 
available to the pion-nucleon two-body subsystem 



E vN = E- \l + k N ) 2 + m\_ x , (4) 

and then defining the invariant center-of-momentum energy for this system 

u 2 a = El N - (K + k N f . (5) 

The mass of the A — 1 system, rriA-i, differs from the mass of the A-body target, rriA, 
by a nucleon mass and a binding energy, rriA = rn^-i + win + This energy must 



then be shifted by a 'mean-spectral energy', E ms , which is a calculated number |24j that 
approximately accounts for the interaction of the intermediate A 33 (or the intermediate 
nucleon or other hadronic resonance) with the residual nucleus. This first-order potential 
produces results for energies near the A33 resonance that are in remarkable agreement with 
the data. This is because 0] the sum of the second-order effects (Pauli exclusion, true 
absorption and correlation corrections) cancel amongst themselves. The use of invariant 
normalizations, invariant phase space and invariant amplitudes produces the phase-space 
factors that are present in Eq. ^|. 

The optical potential so defined is then inserted into the Klein-Gordon equation which 
is solved numerically to produce our model-exact calculation. 



B. Eikonal model 



The eikonal model we propose to examine results from first replacing the optical potential 
of Eq. ^| by a local potential and then solving for the scattering amplitude arising from the 
use of the local potential and the eikonal propagator. In order to include the Coulomb 
interaction, we divide the scattering amplitude as 

F(q) = F pt (q) + F CN (q) (6) 

where q is the momentum transferred to the pion, q = k' n — k n , F pt is the point Coulomb 
scattering amplitude, and Fqn is the scattering amplitude calculated from the sum of the 
strong nuclear potential plus a Coulomb correction, which is the difference between the 
point Coulomb interaction and the Coulomb interaction of a finite charge distribution. The 
eikonal approximation gives for F C n{q) 

/•CO 

F CN (q) = iK / bdbJ (qb)e^T CN (E, b), (7) 
Jo 

where k^ is the momentum of the incident pion in the pion-nucleus center-of-momentum 
frame, b is the impact parameter, Xpt is the point Coulomb phase (in the eikonal approxi- 
mation), and Tcn is the profile function defined below. 

In order to incorporate the distortion and energy shift caused by the Coulomb interaction, 
we follow the prescription given in Ref. [2l| and write the profile function as 



T CN (E, b)=l- exp{ixcN(E, b)} 

= 1 - exp{ixN[E, b(l + EV c {b)/kl)\ + t Xc (b) - i Xp t(b)}, (8) 

where V c is the Coulomb potential of a uniformly charged sphere, \ c is the phase shift caused 
by that potential, and the nuclear phase shift \n is written as 

1 f°° 

XN (E,b) = -— dzU eik [E-V c (r) } r] } (9) 

<£K n J— oo 

where 
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r = (62 + 2 2 ) i/ 2) (10) 

and t/eik is the strong potential obtained as follows. First we construct a local potential U 
from the pion-nucleon two-body scattering amplitude and the target density by 

U (E,r) = -AnZ[f p (0)P P + fpV 2 p p ] - AnN[f n (0)p n + ^V 2 p n ], (11) 

7T 7T 

where Z and A" are the proton number and the neutron number respectively. The scat- 
tering amplitude f p (0) Lfn(O)] is the pion-proton [pion-neutron] scattering amplitude in 
the forward direction in the pion-nucleus center-of-momentum frame. The finite range of 
the pion-nucleon interaction is included through lowest nonvanishing order and produces 
the additional terms in Eq. |TT] which are proportional to the derivatives of the scattering 
amplitude as shown in Appendix B. 

In order to investigate the importance of the lowest-order corrections to the eikonal 
propagator as derived by Wallace fTS], next we replace the potential U a (E, r) by the potential 



Tj2 nip. J 

UME,r) = U + ^(l + --lnU ). (12) 
The details of the derivation of this correction and higher order corrections can be found in 



Ref. [161. 



III. RESULTS 

In this section we investigate the importance of some of the ingredients of our model- 
exact calculation of high-energy pion-nucleus scattering. We will then examine how well 
our simple eikonal theory reproduces the results of the model-exact theory. 

A. Contributions from nonlocalities in the optical potential 

The most difficult part of evaluating the optical potential in the model-exact theory is 
performing the Fermi-averaging integration. This is known to be very important in the 
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region of the A33 resonance, where proper accounting for the recoil and propagation of the 
pion-nucleus resonance requires a careful treatment of the dependence of the energy uj a 
in Eq. on the momentum of the nucleon. Resonance propagation manifests itself as a 
nonlocality in the optical potential, and this represents the first source of nonlocality whose 
importance we want to assess for high-energy pions. 

To investigate this we calculate elastic scattering of ir~ from 40 Ca first with full Fermi 
averaging and then utilizing a closure approximation. The closure approximation that we 
use results from taking 



"a = - K and W 7TN = E- yjkl + m\_ x (13) 

This is equivalent to setting the momentum of the struck nucleon to zero in the pion-nucleus 
cm. frame. 



On the delta resonance, it has been found p6|,|27|,p8| that the closure approximation is 



totally inadequate and even the 'optimally-factorized' approximation |28] is not quantitative. 
At these higher energies the conclusion is not d priori clear. The pion-nucleon amplitude 
is, partial wave by partial wave, rather energy dependent as the individual partial waves are 
resonant. The total amplitude, however, exhibits only two very broad smooth peaks on a 
large background. 

The resulting differential cross sections for pions at 500 MeV scattering from 40 Ca are 
shown in Fig. We see that the correct treatment of the recoil of the two-body pion-nucleon 
system only affects the depth of the minimum but not to a very significant degree. 

The behavior of the pion-nucleon scattering amplitude off-shell is taken from the 



doorway-resonance model of Ref. Amplitudes which are resonant are separable in 

their dependence on the relative momenta, hence maximally nonlocal. This is the second 
source of nonlocality whose importance we want to investigate. In the doorway-resonance 
model, in each angular momentum, spin and isospin channel, 

(k! I tj LI {u a ) I K) = — (k I t JLI (uJ a ) I K ) — r, (14) 
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where k q is the on-shell momentum, i.e. the momentum corresponding to the center-of- 



momentum energy u a . This nonlocality, which is also factored in the delta-hole model, 
was found quite important |28| for energies on and below the A33 resonance. For numerical 
convenience, we take the form factor to be a Gaussian, v(k) = exp(— k 2 //3 2 ), and we vary (3 
from 500 MeV to 4 GeV to investigate the importance of this nonlocality. The differential 
cross sections for 7r~- 40 Ca scattering calculated with this variation in the form-factor range 
are shown in Fig. |||. We see that the dependence on the off-shell range is weak at this 
energy. 

There exists some confusion in the literature concerning the possible values for f3. First, 
if one does not include invariant phase space factors explicitly, they will appear effectively 
in what one might wish to call the form factor. This difference is then one of semantics. 
However, in deciding what might be a reasonable range over which to vary the form factor, 
one should treat the kinematic factors explicitly. In the P-wave channels, an artificial and 
incorrect increase (in coordinate space) in the range of the form factor will also result if the 



pion-nucleon pole term [E9 is not included in the model of the pion-nucleon amplitude. For 



our purposes here, the difference of interest is in examining how the cross section changes in 
going from the physically motivated doorway-model form factors to the higher momentum 
cutoff. The 4 GeV cutoff produces an approximately zero-range interaction in coordinate 
space. Figure § demonstrates that the nonlocalities contained in a finite range interaction 
do not alter significantly the predicted elastic cross sections for these high energy pions. 

We have found that the nonlocalities in the optical potential at the high energies do not 
have an important effect. This suggests that the simple eikonal approach outlined above 
could prove to be an adequate model for quantitative work at these energies. To investigate 
this, we picture the differential cross sections predicted by each model for 7r ± - 12 C and 7r =t — 
40 Ca at 800 MeV/c in Fig. ^ and Fig. |^. The data are from Ref. Q. As stated earlier, we 
construct the target density from the same wave functions as were used in the momentum- 
space calculation and use the same on-shell pion-nucleon scattering amplitude. For 40 Ca, 
both the location of the minima and the magnitude of the cross section at the forward angles 
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are in good agreement. The agreement is, however, not as quantitative for 12 C where the 
minima are slightly shifted. This difference is presumably caused in part by the fact that 
the kinematics of target recoil enter explicitly in the momentum-space formulation but not 
in the eikonal approach. The eikonal model is thus more quantitative for the heavier nuclei. 
We are examining the eikonal model to see if it can be further improved for the light nuclei. 

B. Important features in the eikonal approximation 

In this section we investigate the role of the Coulomb interaction and the first Wallace 
correction, which are both necessary if one wants to use the eikonal model to quantitatively 
approximate the results of the solution of the model-exact theory. 

First, we find that the Coulomb interaction plays an important role in determining the 
depth of the minima in the differential cross section. To see exactly how important Coulomb- 
nuclear interference is, we compare the calculation of the 7r~- 40 Ca differential cross section at 
a series of energies as shown in Fig. |6] with and without the Coulomb interaction. The Wallace 
correction is not included in either the solid or the dashed cureves. We find that the first 
minimum in the differential cross sections becomes very deep when the Coulomb interaction 
is included, falling well below 10 -2 mb/sr in the region of 500 to 600 MeV. Without the 
Coulomb interaction, the deepest minimum is much shallower and occurs at 780 MeV. The 
inclusion of the Coulomb interaction in the theory is necessary if a quantitative comparison 
with the data for 7r~- 40 Ca is to be made. Coulomb-nuclear interference is destructive for 
7r~ and constructive for ir + , which produces a less dramatic effect for 7r + - 40 Ca scattering in 
this energy region. 

Wallace has shown in Ref. [16[ that the semi-classical approximation can be improved 
if higher order corrections (known as 'Wallace" corrections) are included. We examine in 
Fig. the importance of the first Wallace correction, Eq. for 7r~- 40 Ca scattering. We see 
in Fig. [7|, in the absence of the Coulomb interaction, that the Wallace correction is dramatic 
on resonance, T n = 180 MeV, where it fills in the deep minimum. At the higher energies, 
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the correction interferes constructively until about 700 MeV filling in the minima somewhat. 
At 780 MeV, however, the Wallace correction interferes destructively with U and produce 
a much deeper minimum. 

In Fig. ^ we turn on both the Coulomb interaction and the Wallace correction, we find 
a fascinating interplay between the two. Notice that at 400 MeV, for example, there two 
corrections are out of phase and tend to cancel. At 680 MeV, however, they are in phase 
and interfere destructively with U producing the very deep minimum. By 780 MeV they 
are both out of phase with each other but in phase with the lowest order term and thus fill 
in the minima. Both the Coulomb and the Wallace correction are necessary if there is to 
be a deep minima at 680 MeV, and not at the neighboring energies. We note that the data 
of Ref. was taken at a pion lab momentum of 800 MeV/c which is very near T n = 680 
MeV. The depth of the first minimum of the differential cross section for 7r~ scattering over 
the energy range 400 to 800 MeV should prove a very sensitive test of the existence, or lack 
thereof, of higher order corrections in the reaction dynamics. 



IV. COMPARISON TO OTHER WORK 



Other work on high-energy pion-nucleus elastic scattering, which utilizes the full Glauber 



theory, can be found in Refs. fli7y30| - In particular, a comparison between the Glauber theory 
and the optical model on 7r ± - 12 C has been made in Ref. [[L7|] and the two approaches were 
found to be in good agreement at the higher energies. We earlier noted that the optical 
potential model there did not include the full Fermi averaging, and it was not clear to what 
extent the agreement found was caused by the approximate treatment of the non-localities 
in the optical potential. Our work clarifies that. We also note that the Glauber theory 
they used is considerably more complicated than our eikonal model. To our knowledge, the 
Wallace corrections have not been included in any of the numerical tests of the eikonal model 
in this energy region. Hence, the important interplay between the Wallace and Coulomb 
corrections that we find has not been previously noted. We also found that the center- 
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of-mass corrections in the target wave functions are important for C but not for Ca 



which is in agreement with Ref. Our work strongly supports their conclusion that pion 
scattering at the higher energies becomes theoretically and calculationally simpler to treat 
at a quantitative level. 



V. SUMMARY, CONCLUSIONS, AND FUTURE PROSPECTS 

We have compared the scattering from the fully microscopic and nonlocal lowest-order 
optical potential to the scattering from a simple local potential in the eikonal approximation. 
We have investigated the importance of several features of the full model. Within the 
framework of our "model-exact" optical model, these include Fermi averaging and the off- 
shell behavior of the pion-nucleon scattering amplitude. For the eikonal theory, we examine 
the importance of the first Wallace correction and find that this correction is responsible for 
a noticeable improvement of the eikonal theory in comparison to the "model-exact" theory. 
In all cases, the Coulomb-nuclear interference is important and must be included if one is 
to compare with data. 

This semi-classical eikonal theory put forth here is very simple, even simpler than the 
Glauber theory utilized by others. It appears to be quantitatively valid at high energy, at 
least for the first several minima in the differential elastic cross section. This makes possible 
a much simpler reaction theory than has been needed at energies at and below the A33 
resonance. 

At high energies and/or for heavy nuclei, use of the full momentum-space optical model 
becomes very time-consuming and sometimes computationally impossible. In contrast, the 
eikonal theory is relatively simple and fast on the computer. It therefore becomes a matter 
of considerable practical importance to realize that the physics is faithfully reproduced by 
the simple version of the theory out to the position of the second minimum in the differential 
cross section. 

An examination of the pion-nucleon cross section shows the complex interplay of reso- 
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nances. Ultimately, it is interesting to explore how resonances behave in the nucleus and 
how the specific features of the optical potential are suited for nuclear structure studies. In 
this pursuit, it will be very helpful to have a simple version of the theory as developed here. 
We will pursue these questions in subsequent work. 

The momentum space optical model has been extended to treat kaon-nucleus scattering 
32fl . Results show qualitative agreement with previous work ||. Modifications to the kaon- 



nucleon amplitudes in the nuclear medium are needed to eliminate discrepancies between the 
theory and the data. On the quantitative level, results from the momentum space approach 
give larger discrepancies between theory and data than do those from the coordinate space 
calculations used in Ref. ||. Our eikonal model, which is also a coordinate space approach 
and been shown to be a good approximation to the "model-exact" momentum space theory 
for high-energy pion scattering, can also be modified to treat the kaon scattering. Results 
from the eikonal model may reveal the sources of the discrepancies between the momentum 
space and the previous coordinate space calculations. Work on the kaon-nucleus scattering 
utilizing the eikonal model is in progress. 
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APPENDIX A: CENTER OF MASS CORRECTIONS 



The scattering theory developed in Ref. [11] is covariant. However, to implement it 



covariantly would require a covariant model of the nucleus. This would clearly be overkill 
as the nucleus for the cases of interest is not moving relativistically. A Galilean invariant 
model would thus suffice. Contemporary models of the nucleus, however, are not Galilean 
invariant. We use Hartree-Fock models of the nucleus which contain spurious center-of-mass 
motion. The approach we adopt is discussed in detail in Ref. |T2j . The Hartree-Fock results 
are fitted to elastic electron scattering. The spurious center-of-mass motion is there treated 
by requiring that momentum conservation holds for average values. We make the same 
approximation for pion-nucleus scattering. 

Specifically we calculate the target wave functions from the Hartree-Fock model and 
then calculate the Fourier transform of the density. We require the intrinsic density, i.e. 
the density without the center-of-mass motion. If the wave function were a pure harmonic 
oscillator, then the correction could be made exactly. Here, we make the correction as if the 



Hartree-Fock result were exactly a harmonic oscillator. The result [31] for obtaining the 
intrinsic proton density (f cm from the Hartree-Fock density p p HF is 

tiJr) = J d 3 r 1 (B^)- 3 e^p p HF (\r--n\) , (Al) 

where B = b/y/A and b is the usual harmonic oscillator parameter. The Hartree-Fock 
calculation uses a slightly different b for each orbital so we use the average value in making 
the center-of-mass correction. 

We use pP. m to calculate the r.m.s. radius and compare this to the r.m.s. radius measured 
in electron scattering (with the finite size of the proton removed). For 12 C and the Hartree- 
Fock result from Ref. J21] we find that the r.m.s. radius is 2.365 fm from the Hartree-Fock 
calculation while it is measured to be 2.297 fm in electron scattering. We therefore scale the 
&'s in the Hartree-Fock calculation by .9715 (a 3% correction) so that the correct measured 
r.m.s. radius is reproduced. After the scaling, we calculate the predicted electron scattering 
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and find that these wave functions reproduce electron scattering data out to a momentum 
transfer which is larger than is encountered in pion scattering. For 40 Ca, these corrections 
were already made by Negele [^UJ and his wave functions are consistent with the electron 
scattering data. 



APPENDIX B: FORWARD SCATTERING AMPLITUDE 

Using a Taylor expansion, we can write the forward scattering amplitude for small angles 

as 

df(6) 



f(0) « /(0) + 



dcosQ 



(cos 9- 1) 



/(0)"f(0)y, (Bl) 



where 



AO) im 



dcos9 

E/»(°) 



cos 9=1 

d Pi (cos I 



dcosO 



cos 9=1 



£/»(0)^±^. (B2) 



2 



We further use the approximation 



9 « (B3) 

and the fact that the momentum transfer q in momentum space corresponds to zV in coor- 
dinate space, to rewrite Eq. [Bl] as 



f( e ) « /(0) + ^#V 2 . (B4) 

7T 

To construct a local optical potential, we multiply f{6) by the target density p and obtain 
Eq.0. 

The forward scattering amplitude /(0) and its derivative /'(0) are converted from phase 
shifts in the pion-nucleon cm. frame, then transferred to the pion-nucleus cm. frame 
through the relation: 
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7(0) = ^Aat(O), (B5) 

where k and /„-jv(0) denote the pion momentum and forward scattering amplitude in the 
pion-nucleon cm. frame. 
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FIGURES 

FIG. 1. Profile function S(b) for 12 C, Ca and 208 Pb. The arrows indicate the depth to which 
the pion (labeled by its energy) can penetrate. 

FIG. 2. The differential cross section for elastic scattering of 7r~ from 40 Ca at 500 MeV verses 
center-of-mass angle. The solid curve is the result of a calculation in which the full Fermi averaging 



is performed exactly, while the dashed curve uses the closure approximation defined in Eq. 13. 



FIG. 3. The differential cross section calculated from different form factor ranges (5 defined 
in Sec. Ill for n~- Ca elastic scattering at 500 MeV. The solid curve is from the "model-exact" 



theory. The doted curve corresponds to j3 = 4 GeV. The dashed curve is the result from setting 
P = 500 MeV. 

FIG. 4. The elastic differential cross sections for 7r ± - 12 C scattering at a pion lab momentum 
800 MeV/c. The solid curves are from the momentum space calculations, while the dashed curves 
are from the eikonal model. Data are obtained from Ref. Hi. 



FIG. 5. The same as Fig.|| except the target is 40 Ca. 

FIG. 6. 7r - - 40 Ca elastic differential cross section calculated from the eikonal model at labeled 
pion lab energies. The Coulomb interaction is included in all the solid curves, but not the dashed 
curves. The Wallace correction is not included here. 

FIG. 7. 7r~- 40 Ca elastic differential cross sections calculated from the eikonal model at labeled 
pion lab energies. The dashed curves are the same as in Fig. g The solid curves include the 
Wallace correction but not the Coulomb interaction. 

FIG. 8. The dashes curves are the same as in Figs. || and [?| The solid curves now contain both 
the Wallace correction and the Coulomb interaction. 
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